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1 Introduction 



In the present work we discuss phenomenological applications of the formahsm presented in 
ref. [1], relevant to the production of heavy quarks in photon-hadron collisions. All technical 
details and theoretical motivations for the formalism are discussed there. In the formalism of 
ref. [1] one obtains a cross section which is accurate to 0{aemCig) (that is, to the next-to-leading 
order level) and which, in the large Px region, includes consistently all enhanced logarithmic 
terms of the form acmas(as logpT/mY and aemOs (o^s logpT^/my. Terms of order ©(aemOs) with 
no logarithmic enhancement and terms of order cx^^^cx •^{aslogpT/my are not included. Thus 
this approach improves both the fixed order (improperly called massive) and the resummed 
(improperly called massless) calculations. 

Our starting point is the "matched" formula [1] 



where FONLL stands for fixed-order plus next-to-leading logarithms, and 

• FO is the fixed order, 0{aemCts) result; 

• RS is the resummed result, which includes all terms of the form Qemas(«s logPT/^T-)' and 
c^emC^si'^s ^ogp-r/my, and neglects all terms suppressed by powers of the heavy quark mass 
m; 

• FOMO is the massless limit of FO, in the sense that all terms suppressed by powers of m 
are dropped, while logarithms of the mass are retained; thus FOMO is the truncation of 
RS to order aemOs! 

• G{m,pT) is an arbitrary dumping function, that must be regular inpT, and must approach 
1 up to terms suppressed by powers of m/px at large p-j-; our standard choice is 



which for c 7^ dumps the RS — FOMO term in the small-px region. 

We remind the reader that in the description given above the hadronic photon contribution 
is also included; the counting of the orders in ofem and remains the same, provided one 
counts the parton density functions of the photon as carrying a Oem/as factor. A more detailed 
discussion can be found in ref. [1] . 

This paper is organized as follows. In sect. 2 we specify our choices for the input parameters 
and the parton density functions (PDF's). Furthermore, the specific cuts that the Zeus and HI 
collaborations use in their analyses are reported in this section. All the subsequent results will 
be based upon these choices. 

In order to present meaningful phenomenological results, we need to consider a non-perturbative 
fragmentation function (NPFF) to account for the hadronization of the charm quark into a 
charmed hadron. The NPFF is the source of several uncertainties affecting the theoretical pre- 
dictions. First of all, various forms of NPFF have been proposed in the literature (see ref. [2] 
and references therein for an updated review). Furthermore, while the application of NPFF is 
unambiguous in the large p-r region, at moderate and small p-r there are ambiguities, depending 



FONLL = FO + (RS - FOMO) x G(m,pT) 



(1.1) 



G{m,pT) = 




(1.2) 
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upon the definition of the z variable (i.e., whether it refers to energy, momentum or "plus" 
component), and upon the choice of the frame in which it is defined. The impact of these 
ambiguities will be assessed in sect. 3. 

Phenomenological results have also uncertainties depending upon the accuracy of the pertur- 
bative calculations (which we estimate by varying the renormalisation and factorisation scales) , 
the allowed range in the value of the heavy quark mass and Aqcd, and the choice of the PDF's. 
These uncertainties will be discussed in sect. 4. 

The HI [3] and the Zeus [4] collaborations have studied photoproduction of D* mesons. The 
two experiments have different acceptances and apply different cuts to their data samples. In 
sect. 5 we compute the distributions presented in refs. [3, 4], implementing the same cuts that 
are used there. 

In sect. 6 we give our conclusions. 

2 Input parameters 

All the results presented in this paper are relevant to photoproduction experiments where the 
photon arises from collinear radiation off an incoming electron. The radiation spectrum depends 
upon the cuts that are applied to define the photoproduction regime, which are typically 
cuts. Further kinematic cuts are applied either to the the variable y (the ratio of the photon 
energy over the incoming electron energy) or, which is equivalent, to the variable W (the photon- 
proton center-of-mass energy). The cuts used in refs. [3, 4] are summarised in table 1. The HI 





Q^ax (GeV^) 


?/min 


J/max 


Zeus 


1 


0.187 


0.869 


HI ETAG44 


0.009 


0.02 


0.32 


HI ETAG33 


0.01 


0.29 


0.62 



Table 1: Q"^ and y cuts used by HI and Zeus in refs. [3, 4] to define a photoproduction event. 

collaboration uses two different electron taggers, denoted as ETAG33 and ETAG44, and has 
thus two samples of photoproduction data. The electron and proton energies are fixed to 27.5 
and 820 GeV respectively. The y range quoted for Zeus corresponds to 130 <W < 280 GeV. 

We compute the photon spectrum in the Weizsacker- Williams approximation, in the im- 
proved form of ref. [5], which includes non-logarithmic terms enhanced at small y. 

The central values and ranges of the physical parameters that we adopt in the present paper 

(5) 

are summarised in table 2. The range in Aqcq reported in the table corresponds to the range 
a^iM'z) = 0.118 ± 0.002. We use the two-loop expression of throughout this work. The 





low 


central 


high 


Proton PDF 


CTEQ5M [6] 


Photon PDF 


AFG [7] 


Charm quark mass m 


1.2 GeV 


1.5 GeV 


1.8 GcV 


■''■QCD 


0.203 


0.226 


0.254 



Table 2: Central values and ranges of the physical parameters adopted in this paper. 
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factorization and renormalization scales {fip and fin) are set equal to /xo/2, /xq and 2/io, where 
Mo = A/plT+rT? is our default choice. 

In table 3 we report the kinematic ranges of the cross section measurements used in the 
Zeus and HI analyses. Notice that the D* kinematics are described by Pt, rj (pseudorapidity) in 
the Zeus analysis, and hy pT,y (rapidity) in the HI analysis. Unfortunately, the same notation 
(y) is traditionally used to indicate the Weizsacker- Williams variable, and the meson rapidity. 
However, no confusion should possibly arise. In what follows, y will always indicate the meson 
rapidity. 



Zeus 


Pt{D*) > 2 GeV 


\v{D*)\ < 1.5 


HI ETAG44 


PriD*) > 2 GeV 


< 1.5 


HI ETAG33 


Pt{D*) > 2.5 GeV 


\y{D*)\ < 1.5 



Table 3: Visible ranges for the HI and Zeus analyses of refs. [3, 4]. 



3 Non-perturbative fragmentation function 

We begin by discussing the ambiguities related to the implementation of the non-perturbative 
fragmentation in our calculation. A NPFF has to be convoluted with our full FONLL result. 
The convolution is unambiguously defined only at large transverse momenta (i.e. Pt » m). 
One writes 

= / d^k dzDMz)'^^5^ {k-zk) , (3.1) 

where is the cross section for the production of the heavy-flavoured hadron^ H with momen- 
tum k, and ctq is the cross section for the production of the heavy quark Q with momentum k. 
-Dnp(-z) is the NPFF. In eq. (3.1) we assume that fragmentation scales the 3-momentum of the 
heavy quark. Different prescriptions are however possible: for example, one could assume that 
either k^ ov + \k\ are scaled instead of k. In all cases, it is assumed that k remains parallel 
to k, and the mass-shell condition is imposed to fix the remaining components of the meson 
momentum. Only in the large-px limit all these prescriptions coincide. Also, none of these 
definitions is boost invariant, since the condition of parallelism is frame dependent. Again, in 
the large-px limit boost invariance is recovered. 

For the meson momentum k we expect the condition k^ = to be satisfied. However, for 
the purpose of computing distributions, we can also use k^ = m?, and consider the difference 
to be just another ambiguity in our result. The reader may find unreasonable to adopt this 
second choice, since physically the on-shell condition of the meson H must apply. However, in 
the context of a pcrturbative QCD computation, setting k'^ = mjj may appear unreasonable 
from other points of view. For example, if we produce a quark with maximal Pt, and fragment 
it with z = 1, setting k^ = mjj would yield a meson with the same p^ but larger energy, which 
violates energy conservation. Similar examples can be given at small px and large rapidities. If 
we are far from the phase space boundaries, which is always our case, both choices are sensible, 
and their difference can be considered as an inherent uncertainty of the procedure. 

^In what follows, we assume the hadron J? to be a D* meson, since the HI and Zeus data that we use in this 
paper are relevant to D* mesons. 
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In the following, we will always use the Peterson fragmentation function [8] 

-2 



1 



z 1 — z 



(3.2) 



For the sake of the studies performed in the present section we will take it normalized to 1, i.e. 
we do not include the D* production fraction. Although other choices of fragmentation function 
are possible, we only use eq. (3.2) because it is the most commonly adopted, and because recent 
phenomenological fits of the e parameter, adequate to the present work, do exist [9, 10]. In 
refs. [2, 11] the reader will find studies of alternative forms. 

Our default choices for the non-perturbativc fragmentation cirG clS follows: 



the z variable is defined by the equation k = zk; 



the fragmentation frame (i.e., the frame where the equation k = zk holds) is the HERA 
laboratory frame; 



• the mass shell condition is k 



m 



2. 



• the value of the Peterson parameter is e = 0.02. 

We now study the effect of the ambiguities in the definition of the non-perturbative frag- 
mentation on the pseudorapidity and on the transverse momentum distributions for the Zeus 
kinematic cuts. Our conclusions would not change had we used HI cuts. We consider first 
the effect of the change of the fragmentation frame. We call scheme (A) our default choice 
of frame. We call scheme (B) the alternative choice of a longitudinally-boosted fragmentation 
frame, boosted with rapidity equal to the quark rapidity. Thus, in this frame the quark rapidity 
is zero. We notice that in scheme (A) the quark and the meson have the same pseudorapidity, 
while in scheme (B) they have the same rapidity. 
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Figure 1: Fragmentation scheme and e dependence of the pseudorapidity and transverse momentum 
distributions with Zeus cuts. The lower plots show the ratio of the variation over the default choice. 

In figure 1 we show the difference in the pseudorapidity and transverse momentum distri- 
butions when using schemes (A) and (B). Furthermore, the effect of varying the value of the 
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Peterson parameter e (in scheme (A)) is also shown. As we can see, the dependence on the frag- 
mentation scheme is visible only for jj-r < 3 GeV. The rapidity distribution is mostly affected 
because of the relatively low cut. The effect is generally modest, below 5% in the range of 
interest. When increasing e, the spectrum becomes softer. The spectrum with e = 0.036 is 
below the one with e = 0.02 for Pt > 1.5 GeV; in the tail, the difference is of the order of 30%. 
As in the case of the fragmentation frame dependence, the behaviour of the p^ spectrum affects 
the ri distribution through the cut. When changing e from 0.02 to 0.036, the rj spectrum is 
almost unchanged in shape, and the normalisation decreases of about 10%. 
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Figure 2: Effect of the on-mass-shell prescription in the pseudorapidity and transverse momentum 
distributions with Zeus cuts. The lower plots show the ratio of the variation over the default choice. 

In figure 2 we show the difference in the pseudorapidity and transverse momentum distri- 
butions when using the mass of the quark (our default choice) for the on-mass-shell condition, 
as opposed to the mass of the D* meson (mu* = 2.010 GeV). The comparison is performed in 
scheme (B), since in scheme (A) there is no difference in pj- and 77 for the two mass choices. The 
effect is only visible below Pt = 5 GeV. In the case of HERA experiments, only px > 2 GeV 
is relevant; in this region, the variation is below 10%. Furthermore, larger relative variations 
are observed only for very small p^, where the differential cross section is actually smaller in 
absolute value. 

4 Dependence upon input parameters 

We now turn to the study of the sensitivity of our result to the chosen range of parameters. We 
begin by showing in fig. 3 the sensitivity to the choice of Aqcd- As we can see the sensitivity is 
quite small, below 10%. This is a consequence of the fact that the allowed range for Aqcd has 
become fairly narrow at present. 

We do not include in our study the correlation of Aqcd with the parton densities, since PDF 
sets with flavour matching conditions consistent with our scheme that allow for Aqcd variation 
are not available at present. For the same reason, we do not include a study of the uncertainty 
due to the choice of the PDF set. In general, the inclusion of such an uncertainty is a highly 
non trivial problem, and it is far from having a satisfactory solution at present. However, it 
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Figure 3: Sensitivity of the transverse momentum and pseudorapidity distributions to the choice of 



must be said that current PDF sets have rather similar gluon densities in the x region probed 
by charm production at HERA; therefore, we do not expect the PDF dependence to be a major 
source of uncertainty in our case. 
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Figure 4: Sensitivity of the transverse momentum and pseudorapidity distributions to the choice of 
the charm mass. 

The mass sensitivity is shown in fig. 4. The cross section is most sensitive to the mass value 
for Pt of the order of 1 GeV. In the visible range {pt > 2 GeV) the sensitivity is modest, of the 
order of 10%. 

It has to be pointed out that our study of mass sensitivity is somewhat incomplete. In fact, 
we do not have the possibility of varying the charm mass in the evolution of the parton densities. 
Studies performed in the hadroproduction case have indicated that these effects are small. 

The renormalization and factorization scale dependences are displayed in figures 5 and 6 
respectively. The factorization scale choice has little impact on our results, except for small 
transverse momenta. In particular, when fip = /io/2, at small the cross section becomes 
unphysical, since it depends upon the parton densities at low scales, a region where current 
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Figure 5: Sensitivity of the transverse momentum and pseudorapidity distributions to the choice of 

the renormalization scale. 




Pt (GeV) r) r? 



Figure 6: Sensitivity of the transverse momentum and pseudorapidity distributions to the choice of 
the factorization scale. 

parametrizations are unreliable. We simply truncate the theoretical predictions at small trans- 
verse momenta in this case. The renormalization scale dependence is instead sizable in the whole 
production range. Since no data are available for p-j- < 2 GeV, the inspection of figures 3-6 allows 
us to conclude that the renormalization scale dependence is the largest source of uncertainty 
when comparing the theoretical results to the data. 

5 Comparison with data 

In this section, we compare the theoretical predictions, obtained with our FONLL formalism, 
with the data relevant to D* production, presented by the HI and Zeus collaborations in refs. [3, 
4]. For the sake of this comparison, we now need to multiply the NPFF of eq. (3.2) by a factor 
of 23.5 %, to take into account the c D* fraction [12]. Furthermore, we multiply by a 
factor of 2 in order to get the D*~^ + D*~ cross section, which is what the experiments quote 
in refs. [3, 4]. HI data are presented in ref. [3] as 'yp cross sections, obtained by dividing the 
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measured photoproduction ep cross sections by the Weizsacker- Williams flux; here, we quote 
the data as ep cross sections. 

In QCD, it is very often the case that the theoretical uncertainties are larger than the 
experimental errors. Thus, an estimate of these uncertainties is mandatory for a meaningful 
comparison with the data. Usually, one takes the maximum and minimum values of the theo- 
retical predictions, in the allowed range of the input parameters, as an estimate of the upper 
and lower theoretical errors. In our case, since the computation is numerically intensive, this 
procedure is too time consuming. We have thus adopted the following procedure. We construct 
an error band for our predictions by adding linearly all the variations that we have discussed in 
the previous section. In other words, for each parameter x (mass, Aqco, A^r, A*f) that has central 
value xo and upper and lower values Xh and x; respectively, we define an upper and lower error 

(5^P = mayi[a{xo),a{xh),cr{xi)] - a{xo) , (5.1) 
5^-" = aixo)~mm[aixo),aixh),cTixi)], (5.2) 

with all the other parameters in the cross section kept fixed to their central values. We then 
define our upper (lower) end of the band as the central cross section plus (minus) the sum of all 
upper (lower) errors. By interpreting eqs. (5.1) and (5.2) as true errors affecting the theoretical 
predictions, we could have also adopted the strategy of adding them quadratically. We refrain 
from doing this, since it is actually very doubtful that they obey a Gaussian law. However, it 
must be stressed that, in practice, the quadratic sum differs only by a tiny amount from the 
linear one, since they are both dominated by the effect due to the renormalization scale, as 
discussed in the previous section. 




Pt (GeV) n 



Figure 7: Comparison of Zeus data with the theoretical predictions. The dotted curve is the central 
value, the band included between the solid curves is obtained by adding linearly the uncertainties, the 
dashed band does not include the factorization scale uncertainty. 

The lower error due to the factorization scale dependence is undefined for small momenta. 
Thus, the lower end of the band is undefined for low momenta if the factorization scale depen- 
dence is included in the error. We therefore present one band that includes the factorization 

scale error, and one band that does not include it. Since the factorization scale dependence is in 
general rather small, the difference between the two bands is non-negligible only in the small-px 
region, basically outside the region visible to experiments. 
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We begin with Zeus data from ref . [4] , for which we present a comparison with our FONLL 
predictions in fig. 7. As often stated in the Hterature, the data he above the theoretical pre- 
dictions obtained with the default choice of parameters. On the other hand, we see that the 
data are marginally consistent with our upper band. However, the data seem to indicate a 
Pt spectrum harder than that suggested by QCD. As far as the pseudorapidity distribution is 
concerned, QCD appears to do a decent job, except in the positive-ry region. The rightmost 
data points, however, are by far the ones affected by the largest errors. 




Pt (GeV) 



Figure 8: As in fig. 7, for the HI ETAG33 data. 
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Figure 9: As in fig. 7, for the HI ETAG44 data. 

The data from HI are compared to our calculations in fig. 8, for the ETAG33 sample, and 
in fig. 9, for the ETAG44 sample. In the case of the y spectra relevant to the cuts 2.5 < p-r < 
3.5 GeV, 3.5 < px < 5 GeV, and 5 < p-r < 10.5 GeV, we multiply the data by the width of 
the Pt bin size, since they are originally quoted in ref. [3] as da/dydp-r cross sections. The p-r 

spectra arc perfectly described by our FONLL predictions; all the data points are consistent 
with the default predictions within one standard deviation, with the marginal exception of the 
largest-pT point in the ETAG44 sample. Shape-wise, the y spectrum of the ETAG33 sample 
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is also statistically compatible with QCD predictions; the same can be said for the ETAG44 
sample, although in this case the statistical significance is clearly lower. 

In order to assess more quantitatively the size of the agreement (or disagreement) between 
theory and data, wc now present a comparison between the two in a different form. We compute, 
for each data point, the corresponding FONLL default prediction, and we consider the ratio 
data/FONLL. The results are presented as open points in fig. 10 for the spectra, and in 
figs. 11 and 12 for the ij and y spectra. The figures also display bands, whose pattern is 
identical to that of the bands in figs. 7-9. The current bands are in fact obtained by computing 
the ratios of FONLL bands over FONLL default predictions. 




Pt (GeV) 

Figure 10: Open points: ratio of data over FONLL default predictions, for transverse momentum 
spectra measured by Zeus and HI. Solid and dashed curves are the ratios of FONLL bands over FONLL 
default predictions. 
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Figure 11: As in fig. 10, for Zeus r] distribution (left) and HI ETAG33 y distribution (right). 
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Figure 12: As in fig. 10, for HI ETAG44 y distribution. 



The lower inset of fig. 10 clearly shows that four (out of six) data points relevant to the Pt 
spectrum measured by Zeus are more than one standard deviation away from the upper end of 
the theoretical band. The measured spectrum appears to be harder than predicted; however, 
the effect is probably statistically not significant in the case in which the lowest-px measurement 
is excluded from the data sample. All but the lowest-px point are more than 50% away from 
the default prediction, but still within a factor of 2. On the other hand, there appears to be a 
nice agreement between HI ETAG33 and FONLL predictions; ETAG44 data are also reasonably 
reproduced, but there the large errors prevent from reaching firm conclusions. 

The case of the 77 and y spectra is dealt with in figs. 11 and 12. As far as the Zeus data are 
concerned, fig. 11 confirms the trend already visible in fig. 7: for all the px cuts, the measured 
distributions are more enhanced in the positive-77 region than FONLL predicts. This trend is 
however weaker if the rightmost data points, those affected by the largest errors, are excluded 
from the data sample. Also, by excluding these points all data are within a factor of 2 from 
the default theory result. In general, most of the experimental results lie inside the theoretical 
uncertainty band. As already in the case of the spectrum, HI data do not suggest any 
statistically-significant deviation from FONLL predictions. 

From this study, we see that Zeus and HI data compare differently to FONLL predictions, 
the former displaying some disagreements with theory that are not present in the latter. It is 
difficult to understand whether the two data sets are compatible within errors, since the two 
experiments have different visible regions, and use different observables (Zeus use rj, HI use y). 




Figure 13: Comparison of Zeus data with the theoretical predictions. The solid curve is FONLL, the 
dashed curve is FO (at NLO). 
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The error bands presented in this section do not include any of the uncertainties relevant 
to the non-perturbative fragmentation discussed in sect. 3. This appears to be justified for the 
choice of the fragmentation frame and of the on-mass-shell condition, given the size of the effects 
induced. The choice of the e parameter is phenomenologically more relevant, since it can give 
an up to 20% effect on the Pt spectrum. However, it has to be kept in mind that e is not a 
physical quantity; in particular, its value crucially depends upon the scheme in which the hard 
cross section, which is convoluted with the NPFF, is defined. According to ref. [10], e = 0.02 
is the value to be used in the case of a FONLL computation in order to obtain a meaningful 
comparison with data. The study of sect. 3 simply suggests that a larger (smaller) e value gives a 
softer (harder) spectrum. Thus, the comparison of FONLL predictions to Zeus data appears 
to imply the necessity of an even smaller e value than the LEP data fit performed in ref. [10] 
suggests. On the other hand, HI data appear to be consistent with e = 0.02. 

It is interesting to compare the FONLL and FO calculations. The two predictions are 
displayed in fig. 13 for the Zeus cuts. FONLL results are obtained with the default parameters. 
We show two sets of FO curves. One set (dotted curves) is obtained with default parameters 
except for the value of e, set equal to 0.036. This is in fact the value that gives the best fit to 
e+e~ D* production data when a FO calculation is used [10]. The other set (dashed curves) has 
e = 0.02. We include it in order to emphasise, in comparison to the FONLL result, the effect 
of the large-log resummation, which is only present in the latter. As expected (see sect. 1 and 
ref. [1]), FONLL and FO predictions (with the same e) coincide at small Pt^s; the difference 
between the two becomes visible for p-r > 4 GeV. Where the two predictions differ, FONLL must 
be considered superior. Notice that although none of the two predictions agree well with data, 
the shape is better described by the FONLL calculation. We interpret this result as evidence of 
the onset of resummation effects in the data. In fact, the FO result requires larger value of e in 
order to compensate for the lack of resummation effects. 

We do not present here RS predictions, 
since we believe that in most of the visible 
range the transverse momentum is simply 
too small for them to apply. The RS re- 
sult is below our FONLL result by 20% to 
30% in the region 4 < Pt < 20 GeV, the 
difference becoming smaller as p^ increases. 
This difference is due to our choice of the c 
parameter that appears in eq. (1.2), which 
we fixed to 5 on the basis of the fact that, 
in the hadronic component of the cross sec- 
tion at the NLO level, the massless limit of 
the cross section is a good approximation of 
the massive cross section only for ^ 5m 
[1, 13]. We show our prediction for the ex- 
treme choice c = in fig. 14. The c = choice has a softer spectrum, and it undershoots our 
default choice by at most 25% in the large region. 

The definition of a photoproduction event in Zeus and HI is slightly different, since the 
allowed photon virtuality is much smaller in HI. Thus, one is tempted to ask whether the 
differences of HI and Zeus data compared to QCD calculations are due to the fact that Zeus 
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is not completely in the photoproduction regime. The Weizsacker- Williams approximation is 
valid up to non-factorizable terms of order [5] 
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which yields 5 ~ 0.1 for Zeus, and 5 ~ 0.001 for HI. This does not mean that one can expect 
non-factorizable terms of order 10% in the Zeus case, since the leading Weizsacker- Williams term 
is logarithmically enhanced by a factor of log4m^/mg ~ 17. Thus, one expects non-factorizable 
terms to be of the order of a percent in the Zeus case, and one hundred times smaller in the HI 
case. It is thus unlikely that they could be the cause of the differences. A direct comparison of 
Zeus and HI data, using the same kinematical cuts, would thus be useful to understand whether 
the two data sets are consistent among each other. 



6 Conclusions 

In this paper, wc have presented comparisons between theoretical FONLL predictions and data 
relevant to the photoproduction of D* mesons, as measured at HERA by the Zeus and HI 
collaborations. We have studied the uncertainties affecting QCD predictions, and thus estimated 
the largest possible range for theoretical cross sections. Generally speaking, the agreement 
between theory and data is acceptable, after a tuning of the input parameters. While HI data 
do not display any significant discrepancies with FONLL predictions, Zeus data seem to suggest 
a harder transverse momentum spectrum, and an excess towards the positive-pseudorapidity 
region. A detailed comparison between the two sets of data is not possible due to the different 
observables and visible regions used by the two experiments. 

Having established that QCD does a reasonable job in describing charm data, it is inter- 
esting to notice that Zeus and HI collaborations have now presented several results on bottom 
production [14, 15, 16, 17], that compare much worse to QCD predictions than analogous charm 
results. Since QCD expectations for charm agree reasonably with data, one would expect even 
better agreement for the bottom quark, because of the larger mass. It is thus unlikely that the 
bottom excess may be attributed to a failure of QCD predictions. 
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